Sensors 2009, 9, 5952-5967; doi:10.3390/s90805952 



OPEN ACCESS 



sensors 

ISSN 1424-8220 

www.mdpi.com/journal/sensors 

Article 

Oscillation Control Algorithms for Resonant Sensors with 
Applications to Vibratory Gyroscopes 

Sungsu Park '*, Chin- Woo Tan 2 , Haedong Kim 1 and Sung Kyung Hong 1 

1 Department of Aerospace Engineering, Sejong University, Seoul 143-747, Korea 

2 PINC Solutions, Berkeley, CA 94709, USA 

* Author to whom correspondence should be addressed; E-Mail: sungsu@sejong.ac.kr; 
Tel.: +82-2-3408-3769; Fax: +82-2-3408-4333 

Received: 30 June 2009; in revised form: 22 July 2009 / Accepted: 27 July 2009 / 
Published: 27 July 2009 

Abstract: We present two oscillation control algorithms for resonant sensors such as 
vibratory gyroscopes. One control algorithm tracks the resonant frequency of the resonator 
and the other algorithm tunes it to the specified resonant frequency by altering the 
resonator dynamics. Both algorithms maintain the specified amplitude of oscillations. The 
stability of each of the control systems is analyzed using the averaging method, and 
quantitative guidelines are given for selecting the control gains needed to achieve stability. 
The effects of displacement measurement noise on the accuracy of tracking and estimation 
of the resonant frequency are also analyzed. The proposed control algorithms are applied 
to two important problems in a vibratory gyroscope. The first is the leading-following 
resonator problem in the drive axis of MEMS dual-mass vibratory gyroscope where there 
is no mechanical linkage between the two proof-masses and the second is the on-line 
modal frequency matching problem in a general vibratory gyroscope. Simulation results 
demonstrate that the proposed control algorithms are effective. They ensure the proof- 
masses to oscillate in an anti-phase manner with the same resonant frequency and 
oscillation amplitude in a dual-mass gyroscope, and two modal frequencies to match in a 
general vibratory gyroscope. 

Keywords: oscillation control; resonant sensor; vibratory gyroscope; averaging method; 
automatic gain control; phase-locked loop 
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1. Introduction 

Robust operation of a resonant sensor, such as a vibratory gyroscope, requires the resonator to be 
driven at resonance with constant amplitude. However, both the amplitude and resonant frequency can 
vary due to environmental factors, such as changes in temperature or stiffness aging. Therefore, some 
form of oscillation control is needed to track the constantly changing resonant frequency or to tune it 
to a specified frequency while keeping the amplitude constant. 

Automatic gain control (AGC) has generally been used to excite the resonator to track the reference 
amplitude. The application of AGC to the drive axis of a vibratory gyroscope is reported in [1], where 
velocity measurement is used to control the velocity amplitude to maintain the reference amplitude. It 
also presented stability analysis results and AGC performance, but did not provide the noise analysis. 
Typically the velocity sensing circuitry produces larger noise than the displacement sensing circuitry 
does. A phase-locked loop (PLL) has been also used to track a resonant frequency in [2], where 
stability and resolution analysis results of PLL frequency tracking system for MEMS fatigue testing 
are reported. However, it did not suggest a method to sustain the specified amplitude. 

There have also been a few studies on controllers which, instead of tracking the resonant frequency, 
tune it to a specified frequency chosen by the designer. The advantages of this method are that it can 
maintain consistent performance, since the sensor can retain the dynamic characteristics regardless of 
environmental factors such as temperature changes, and simplify the signal processing loop that uses 
the resonant frequency as its carrier frequency. In the literature, only Lyapunov based adaptive control 
schemes for vibratory gyroscopes are reported to place the resonant frequency at a specified 
frequency [3,4]. 

This paper presents two algorithms for controlling the frequency and amplitude of oscillation. One 
control algorithm tracks the resonant frequency and the other algorithm tunes it to the specified 
resonant frequency by altering the resonator dynamics. Both algorithms maintain the specified 
amplitude of oscillations. In the first algorithm, AGC and PLL structures are used to control the 
amplitude and track the resonant frequency. The displacement measurement is used to avoid using 
noisy velocity measurement. The second algorithm is similar to [3,4], in that the reference frequency is 
chosen and then resonant frequency is adapted to it. However, it is different in some ways. First, our 
algorithm modifies the PLL structure to tune the resonant frequency and uses conventional AGC to 
regulate the amplitude, making it simpler than an adaptive control approach. Second, our algorithm 
can adapt only the frequency to the reference and set amplitude free rather than regulating it, which 
makes it applicable to mode tuning on the sense axis of a vibratory gyroscope, whereas an adaptive 
control approach cannot. 

The averaging method is used to analyze the stability of the entire feedback system and the effects 
of the displacement measurement noise on tracking and estimation of the resonant frequency. The 
proposed control algorithms are applied to two important problems in a vibratory gyroscope to 
evaluate their performance. The first one is the leading-following resonator problem in the drive axis 
of a MEMS dual-mass vibratory gyroscope without a mechanical linkage between the two proof- 
masses. Adopting the proposed control algorithms enables the implementation of the following 
resonator which precisely traces the oscillation pattern of the leading resonator. The second one is the 
on-line modal frequency matching problem in a general vibratory gyroscope. Applying the first 
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algorithm on the drive axis and the second algorithm on the sense axis enables the resonant 
frequencies of the drive and sense axes to be matched precisely, which significantly improves the 
gyroscope performance. 

2. Frequency Tracking and Amplitude Control 

The equation of motion of a resonator is modeled as a spring-mass-damper system and described as 
a second-order differential equation as shown in (1): 

x + dx + co 2 x = f (1) 

where x is the displacement of resonator, d is the normalized damping coefficient, co„ is the natural 
frequency, and /is control input. 



2.1. Control Algorithm 



The control input used to excite the resonator at constant amplitude, X 0 , while tracking the resonant 
frequency of the resonator is given by: 

/ = ^cos# (2) 

where A is the magnitude of the input and 6 is the instantaneous phase of the input. The magnitude 
should be chosen to make the difference between the measured and the specified amplitude, X 0 , equal 
to zero. Therefore, A should have an AGC form. On the other hand, cos# should have a PLL form 
because it should recognize the resonance when the difference in phase angle between a control input 
and a displacement is 90°. 

The AGC is composed of a rectifier, a low-pass filter, and a comparator [1]. However, in order to 
precisely maintain the specified amplitude, X 0 , we employ a proportional-integral control in the AGC. 
The proposed AGC in this paper is as follows: 

A = K p (X 0 -r) + B 
B = Kj(X 0 -r) 



r = X, 



— \x + n p 



2 



where K p is the proportional control gain; Kj is the integral control gain; r is an estimation of the 
amplitude of x, which is the low-pass filtered value of the absolute value of x; and X a is the corner 
frequency of the low-pass filter. n p is the displacement measurement noise and is assumed to be a 
white Gaussian noise with zero mean and power spectral density (PSD) of <j 2 [m 2 1 Hz] . 

A PLL is composed of a phase detector, a PLL controller, and a voltage-controlled oscillator (VCO). 
Each can be written as follows [2]: 

co = 0 = co Q + K v z 

z = K z y (4) 
y = Z p ((x + n p )cos0-y) 



Sensors 2009, 9 



5955 



where cos# is the output of the VCO. The instantaneous phase angle, 6, is the integral of the 
instantaneous frequency co, which is the sum of a proportional value of the control voltage of VCO, z, 
and the free oscillation frequency, a>o. The VCO control voltage, z, is the integral of the output of the 
phase detector y, which is the calculated difference in phase between the output of the VCO and the 
displacement, and is given as the low-pass filtered value of cos# multiplied by the displacement. In (4), 
K v and K z are the gains of VCO and integral control, respectively, and X p is the corner frequency of the 
low-pass filter. Figure 1 shows a block diagram of the frequency tracking and amplitude control 
composed of an AGC and a PLL. 



Figure 1. Block diagram of frequency tracking and amplitude control. 
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2.2. Stability Analysis 



This section follows the formulation in [1,2] and employs the averaging method to analyze the 
stability of the feedback system depicted in Figure 1 with respect to the design parameters: the control 
gains K p and Ki of AGC, the control gains K v , K z of PLL, and corner frequencies X a , k p of the 
low-pass filter. 

Transforming of the displacement and velocity, x and x , to the amplitude, a, and the phase angle, 
<j> , and applying them to the feedback system composed of (1), and controllers (2), (3), and (4) yield: 



1 



a = 



co 0 +K v z 



[(K p (X 0 -r) + B) sin(0 + </>) cos 6 

+ a{K v K z y + d(a) o +K v z)}sm 2 (0 + <fi) 

+ a{(» 0 + K v z) 2 - co\ }sin(0 + 0)cos(0 + $)] 



to 0 + K v z 



1 



a 



(K p (X 0 -r) + B) cos(0 + <j>) cos 0 



+ {K v K z y + d(co 0 + K v z)}sm(6 + <fi) cos(6 + <f)) 
+ {(co o +K v z) 2 -a 2 „}cos 2 (e + 0)] 



B = K 1 (X 0 -r) 



71 



a cos(# + <f>) 1 —r 



1 n i 



(5) 



(6) 



(7) 
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* = K z y 

y = A {a cos(0 + (jf) cos 0 - y] + An cos 0 



(8) 



In (7), r denotes the estimate of the amplitude of x, so r > 0, and the noise term can be separated by 
using | x + n p | < | x | + | n p |. The mean value of | n p |, which is denoted as no, is calculated as: 



2(7, 



If m„ is defined as: 



n 0 = 



m=\n, 



In 



the mean value of m p is zero, and the PSD of m p and the correlation of m p and n p are given by: 



m'p =n] -n\ ~ 0.36cr; 



m p n p =cr p 



(9) 



(10) 



(11) 



where the bar denotes the stochastic expectation. With (10), (7) can be rewritten 



as: 



r = X. 



^(\acos(0 + <fi)\+n Q )-r 



1 n 
+ A n —m n 

a 2 p 



(12) 



Equations (5)-(8) are simplified to a general stochastic nonlinear state space equation as follows. 



*l = f l( x l) + glW 



(13) 



where: 



x 1= [a <j) B r z yf , w = [m p n p ] 2 
"0 0 0 An/2 0 0 



8i 



0 0 0 0 0 A cosO 



(14) 



A first-order approximation of the Taylor series expansion of the nonlinear function fi(xi) in (13) 
about the mean, x t , yields a time update of expectation as follows: 



x, =f 1 (x 1 ) 



(15) 



Since 9 changes much faster than other variables, the averaging method can be applied to (15), and 
the averaged dynamic equations are obtained as follows: 



X lav _ flav( X lav) 



(16) 



where subscript, av, denotes the averaged value. Equation (16) can be written in detail as follows. 



co 0 +K v z 



±(K p (X 0 -r av ) + B av )sm<t> a 



+ ^ L {K v K z y av +d(a> 0+ K v z av )} 



Sensors 2009, 9 



5957 



K = 



2a„ 



<KJX 0 -r m ) + B„)cos*. 



+ ±{(co 0+ K v z av ) 2 -co 2 n } 



r„„ = X, 



^ av ^ z~y av 

fl_ - _ 
y av =^ P -a av cos</> av -y a 
\ z J 



An equilibrium point of (16) is: 



r 0 = X 0 , a Q =X 0 --n 0 , </) a = -—, B 0 = a 0 dco n , z 0 = — , y 0 = 0 

11 K v 

The Jacobian matrix of the averaged system (16) at the equilibrium point is: 



dx 1 



d 


0 


1 




a 0 dK v 


a 0 K v K z 


2 




2co n 


2co n 


2co n 


0 
0 


d 

2 
0 


0 
0 


0 


~K V 
0 


0 
0 


K 


0 


0 


-K 


0 


0 


0 


0 


0 


0 


0 


K z 


0 


A p a 0 
2 


0 


0 


0 


-K 



The conditions for the above Jacobian matrix to be stable are 

+ A 



d + ^- 



v 



KAfd n \ K, Ad ^ 

2 " 



CO 



n J 



> — — , d 



> a 0 K v K z 



(17) 



(18) 



Equation (18) is the stability criteria of the control system and can be used for selecting control 
parameters for achieving stability. The left-hand side of (18) consists of the design parameters related 
only to AGC, while the right-hand side consists of the design parameters related only to PLL. 
Therefore, we can design each controller separately. If the control system is stable, 0 = co n is achieved 

and the excitation frequency tracks the resonant frequency, co„. In addition, since ao = X 0 - non/2, the 
amplitude converges to the spcified value, X 0 , with a small deviation due to the displacement 
measurement noise. 
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2.3. Resolution Analysis 



As seen in the previous section, the displacement measurement noise not only causes errors in the 
amplitude control, but also affects the resolution of the resonant frequency tracking. Applying the 
covariance propagation equation, this section describes the effects of displacement measurement noise 
on the resolution of the resonant frequency tracking controller. The covariance propagation equation of 
(13) is defined as: 



d 
dt 



(*! -XjXx, -X!) 



(19) 



Expanding above equation using the Taylor series at the mean, xi, and obtaining a first-order 
approximation yields the covariance propagation equation as follows [3]: 



P = P 
i ' i 





T 






+ 








v 5xJ 



^+g,(x)S gl (x) 3 



(20) 



where the PSD of measurement noise vector, 5*, is given by: 

V.JDO" 

s = 



0.36^ a p 



(21) 



Applying the averaging method to (20) yields the covariance equations for (13) at steady state as: 



0 = P l F 1 T +F l P 1 +Q l 



(22) 



where Fi is defined in (17) and: 



a =diag{0, 0, 0, (A a 7r/2Y(036a z p ), 0, ^12} 

The standard deviation, which is the averaged resolution of excitation frequency in the frequency 
tracking control, is derived as: 

^=^^(5,5) (23) 
where P\{5,5) denote (5,5) terms of the numerical solution of (22). 



3. Frequency Tuning and Amplitude Control 



3.1. Control Algorithm 



The frequency tuning control differs from the frequency tracking control in that it adjusts the 
dynamic characteristics of the resonator to match to a resonant frequency specified by the designer, 
instead of its own resonant frequency. This requires controlling the resonant frequency through the 
displacement feedback. Adopting the specified resonant frequency, co s , Equation (1) is rewritten 
as follows: 

x + dx + (o 2 s x + Aco x = f (24) 
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where Aco = co] - co] is the difference between the actual and the specified resonant frequencies. The 

objective of the frequency tuning control is to compensate this value. Therefore, the control input is 
consisted of a part that excites the system at the specified resonant frequency, co s , and another part that 
compensates the difference in the resonant frequency through the displacement feedback. This is 
given by: 

f = Acos(a> f t) + Aa>(x + n p ) (25) 

where A is in AGC form, identical to (3), because it is an amplitude control that maintains the 
amplitude of the displacement at the specified value. Here Ad) is an estimate of Aa>, and n p is the 
displacement measurement noise. The adaptation law of Ad) is derived from the fact that the phase 
difference between the control input and the displacement is -90° when Ad> = Aco . The proposed 
adaptation law in this paper comprises the phase detector and the integral controller and is given by: 



Ad) = K m q 

q = A q ((x + n p )cos(a) s t)-q) 



(26) 



where q is the low-pass filtered signal of the control input multiplied by the displacement, K m is the 
integral gain, and X q is the corner frequency of the low-pass filter. Figure 2 illustrates the block 
diagram of the frequency tuning and amplitude control. 



Figure 2. Block diagram of frequency tuning and amplitude control. 
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3.2. Stability Analysis 



Similarly to the previous section, transforming of the displacement and velocity, jc and x , to the 
amplitude, a, and the phase angle, tf> , and applying them to the feedback system composed of (1), and 

controllers (3), (25), and (26) yields. 



— \(K (X 0 -r) + B) sin(» s * + <j)) cos(<y s 0 

+ Adm p sm(6) s t + (/)) + adco s sin 2 (coj + <fi) 
- aAd>sm(co s t + (j)) cos(co s t + </>)] 



(27) 
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(/> = - 



oo. 



1 



a 



( K P ( X o ~r) + B) cos(oo s t + </>) cos((o s t) 



+ — Abbn cos(oo s t + (/)) + doo s sin(<y^ + <f>)cos(oo s t + </)) 
a 



(28) 



• A<»cos 2 (<y s Z L + ^)] 



B = K I (X 0 -r) 



r*2 



71 



(| a cos{oo s t + </)) | +n 0 )-r 



Abo = -K a q 

q = X q {acos(oo s t + </>)cos(oo s t) - q} + X q n p cos(oo s t) 



(29) 



(30) 



where (29) is identical to (7), and Abo = Aco-Aob denotes the estimation error. Equations (27) (30) are 
simplified to a general stochastic nonlinear equation as: 

x 2 = f 2 (x 2 ) + g 2 w 



(31) 



where: 



x 2 =[a (j) B r Abo q] T , w = [m p n p ] T 

0 0 0 A a x/2 0 0 

A<ysin(»/ + 0) Aobcos(oo s t + <fi)/ a 0 0 0 A cos(oo s t + 0) 



g 2 = 



(32) 



Similarly to the previous section, an approximated time update of expectation of (31) is obtained 
as follows: 



X 2 — ^2 ( X 2) 

Applying the same method yields the nonlinear averaged dynamics of (33) as: 

X 2»v = ^2av( X 2av) 

or in detailed forms as: 



(33) 



(34) 



^-[(K P (X 0 -r av ) + B av ) sin 0 av +doo s a av ] 
loo. 



d> =- 

r av 



loo. 



-(K p (X 0 -rJ + B av )cos<{, av -Aa> av 



B av =K I (X 0 -r av ) 
( 



av a 



— K 

2 



a av + — n 0 ~r av 



1 



a = A —a cosd —a 

"av ? o av 'av ^ a 

V 1 
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An equilibrium point of (34) is: 



r 0 =X 0 , a 0 = X 0 -^n 0 , fa =-y, B 0 = a 0 dco s , A5 0 =0,^=0 
The Jacobian matrix of the averaged system (34) at the equilibrium point is: 



dx. 



d 

2 


0 


1 

2co s 


2co s 


0 


0 


0 


d 
2 


0 


0 


1 

2co s 


0 


0 


0 


0 




0 


0 


K 


0 


0 


-K 


0 


0 


0 


0 


0 


0 


0 




0 


2 


0 


0 


0 


-K 



(35) 



The conditions for the above Jacobian matrix to be stable are: 



r K ^ 
p 



d + 



KAfd K, (d „\ d n K„ 



co 



J 



>—*-, d 
CO, 



2 q 



> 



2co r 



(36) 



If the feedback system is stable, actual resonant frequency is tuned to the specified resonant 
frequency since Acd 0 = 0 is achieved. 

3.3. Resolution Analysis 

Similarly to the previous section, applying the averaging method yields the following steady state 
covariance equation for (31): 



0 = P 2 F 2 + F 2 P 2 + Q 2 



(37) 



where F 2 is defined in (35) and: 



02 = 



2 
0 
0 
0 
0 

Acoa 2 p X q 



2 
0 

0 

0 

0 



0 
0 
0 
0 
0 



0 
0 



0 
0 



V z J 



(0.36^) 0 
0 0 
0 0 



Acoa p X q 



0 
0 
0 
0 



The standard deviation, which is the averaged resolution of the frequency compensation in the 
frequency tuning control, is derived as: 

a 2 = VA(5,5) (38) 



where ^(5,5) denote (5,5) terms of the numerical solution of (37). 
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4. Application to the Drive Axis of a Dual-Mass Gyroscope 

In general, a dual-mass gyroscope has two proof-masses linked to each other by a mechanical beam 
and is designed to oscillate in anti-phase with the same resonant frequency and amplitude [5]. The 
Coriolis force acts along opposite directions for each of the two masses, and an external force acts in 
the same direction, enabling cancellation of the acceleration and other common mode effects. However, 
because of issues such as the difficulties in the drive axis alignment due to the connecting beam, 
possibility of unexpected vibration mode, and high production cost induced by high precision process, 
it has been proposed to remove the connecting beam and have the controller oscillate the two masses 
in anti-phase with the same resonant frequency and amplitude [6]. 

To evaluate the control performance, the proposed control algorithms are applied to the drive axis 
of a dual-mass gyroscope without a mechanical linkage between the proof-masses. The equation of 
motion for the drive axis of a dual-mass gyroscope is described by two second-order differential 
equations as follows: 

x x + d x x x + co x x x = f x 

7 . 2 (39) 

X^ H - it ^^2*^2 — J 2 

where subscript 1 and 2 denote the first and the second proof-mass, respectively, x\, x 2 are the 
displacements of each proof-mass, d\, di are normalized damping coefficients, co\, a>2 are natural 
frequencies, and f\,fi are control inputs. 

The first proof-mass is set as the leading resonator, and the second one is set as the following 
resonator, which precisely traces the oscillation pattern of the leading resonator. The proposed 
frequency and amplitude controls are applied to each mass. The dual-mass gyroscope parameters are 
taken from a prototype fabricated at Sejong University, and the parameters are: 

m x = 2.30 KHz, co 2 = 2.53 KHz 

d x =72.26 sec' 1 , d 2 =79.5 sec" 1 (40) 

Considering likely manufacturing errors, 10% error is added to <x>\ for a>2 calculation, and 10% error 
is added to d\ for J 2 calculation. The PSD of the displacement measurement noise is assumed to be 
<r 2 p = 6.92 x 10 24 m 1 1 Hz . The control law applied to the first proof-mass is frequency tracking and 
amplitude control and is written as: 

f x =A x co$6 x (41) 

where A\ and 9\ are calculated from (3) and (4), respectively. The control law for the second proof- 
mass is frequency tuning and amplitude control. In anti-phase, the second proof-mass should oscillate 
with the same resonant frequency as the first proof-mass, so the control law (25) is modified as: 

f 2 = -A 2 cos 0 X + Acb x 2 (42) 
where 0\ is identical to (41) and Ai and Ad) are calculated from (3) and (26). 

The control parameters are selected to meet the stability conditions in (18) and (36). Table 1 shows 
the values of the control parameters used in the simulations. Note that the values in Table 1 are non- 
dimensionalized based on length 1 um and time \lco\ (sec). 
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Table 1. Non-dimensional values of the control parameters. 



Parameter 


Value 




Xq, K p , Kj, X a 


5,0.018, 1.5 x 


10 5 , 0.05 


a>o, K v , K z , k p 


0.9,0.5, 1.5 x 


10~ 4 , 0.5 




1.5 x 10~ 4 , 0.5 





Figure 3(a) shows that the frequency tracking control on the first proof-mass tracks the resonant 
frequency. The excitation frequency starting from free oscillation frequency of PLL, coo = 0.9&»i, 
reaches the resonant frequency of the proof-mass in 0.4 second and then maintain its value. The 
averaged resolution of the excitation frequency from (23) is ci ~ 0.004 Hz. Figure 3(b) shows the 
response of the estimation error, Aco, showing that it converges to zero around 0.6 second. 

Figure 3. (a) Excitation frequency, (b) Estimation error of the frequency difference. 
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As soon as the estimation error becomes zero, the dynamic characteristic of the second proof-mass 
begins to have the same resonant frequency as the first proof-mass, as observed in Figures 4. The 
estimation error of Aco calculated from (38), is less than 02 ~ 0.01% x Aco . Figure 4(a) is the 
displacement response of the proof-mass. Its amplitude is maintained at the specified value, X 0 = 5 um, 
after some time. Figure 4(b) zooms in the response plotted in Figure 4(a) for 1 msec, showing that the 
two proof-masses oscillate in anti-phase with same resonant frequency and amplitude. 
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Figure 4. Time response of drive axes. 
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5. Application to Mode Matching for a Vibratory Gyroscope 

Most vibratory gyroscopes rely on matching the resonant frequencies of drive and sense axes for 
high performance. However, manufacturing imperfections result in deviations of the resonant 
frequencies from their design values. Therefore, various tuning methods have been developed [7,8] 
and the most commonly employed method is to alter the resonator stiffness by applying bias voltages 
with dedicated electrodes. In this method, the bias voltage should be maintained under the change in 
the operating conditions, such as temperature variations. 

An alternative mode matching method is presented in this section using the proposed oscillation 
control algorithms. This method is adaptive to changes in the environment such that if the resonant 
frequency of sense axis is deviated from that of drive axis, the controller will continually compensate 
the frequency deviation. Therefore it can be used for on-line implementation. 

In this method, the frequency tracking and amplitude control is used in the drive axis, and the 
frequency tuning control in the sense axis. The amplitude control is not used in the sense axis because 
the amplitude should be allowed to change according to the input angular rate. Since the drive axis 
control is same as that in previous section, only the sense axis is considered here. The sense axis of a 
vibratory gyroscope is modeled as: 

y + d y y + co\y = -2Qx + f y (43) 

where y is the sense axis displacement, x is the drive axis velocity, d y is the normalized damping 
coefficient, co y is the natural frequency of sense axis, Q. is the input angular rate about z-axis, and 7^ is 
the control input for sense axis. 

Because the frequency tuning control relies on the displacement measurement, it is required to 
excite the sense axis continually regardless of the presence of an angular rate input. Therefore we 
propose to insert fictitious angular rate into the sense axis, and the control law (25) is modified as 
follows: 
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f y = 2Q 0 X 0 Q) X cos0 x +Aa> y y 



(44) 



where co x and 6 X are drive axis resonant frequency and its integral resulting from drive axis control, 
respectively. Qo is a fictitious angular rate, Xq is the amplitude of drive axis oscillation, and Acb y is the 

estimate of (co 2 Y ~ 0} x)> which is calculated from (26). The first term in (44) is equal to 2Q 0 x if the 

drive axis is at resonance using the drive axis control (41). By substituting (44) into (43), and 
modifying (34) and calculating the averaged equilibrium point yield: 



n -j n 
r o =a 0 +-n 0 , <p 0 =- — , 

2X 0 (Q-Q 0 ) A - 



d. 



, Aco Q =0, q 0 = 0 



where ao is the magnitude of sense axis oscillation, which is proportional to input angular rate. 

Simulations are conducted with the same gyroscope data given by (40), where subscript 1 and 2 are 
replaced by x and y. Fictitious input angular rate is assumed to be 10 deg/sec. Figure 5(a) shows the 
time response of the frequency tuning estimation error. Figure 5(b) shows the time response of the 
sense axis to a step input angular rate of 5 deg/sec applied at 1.2 sec. It is observed that as soon as the 
tuning estimation error becomes zero around 0.8 sec, the dynamic characteristic of the sense axis 
begins to have the same resonant frequency as the drive axis, and the gyroscope is ready to work, with 
the response of sense axis after 0.8 sec larger than that before. Figure 6 shows the estimates of the 
angular rate response to step and sinusoidal input angular rates when the modes are matched. These 
simulation results illustrate that the resonant frequencies of the drive and sense axes are precisely 
matched, and the gyroscope performance is greatly improved. 

Figure 5. (a) Frequency tuning estimation error, (b) Time response of sense axis. 
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Figure 6. Time response of angular rate estimates to (a) step input, (b) sinusoid input. 
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6. Conclusions 

In this paper, two frequency and amplitude control algorithms are presented. One control algorithm 
excites the resonator at its own resonant frequency, and the other alters the resonator dynamics to place 
the resonant frequency at a specified frequency which is chosen by the designer. These control 
algorithms maintain specified amplitude of oscillations. The stability of the entire feedback system was 
analyzed using the averaging method, and the stability criteria were proposed so that it can be used as 
guidelines for selecting the control parameters. In addition, the effects of displacement measurement 
noise on the tracking and estimation of the resonant frequency were analyzed. 

In order to evaluate the performance of the proposed control algorithms, we apply them to two 
important applications. The first one was application to the drive axis control loop for a dual-mass 
gyroscope without a mechanical linkage between two proof-masses. The second one was application to 
modal frequency matching in a vibratory gyroscope for high performance operation. 

Simulation results agreed well with analytical analysis, and demonstrated the effectiveness of the 
proposed control algorithms. In addition, both the analytical and simulation studies showed that it is 
possible to make the two proof-masses of a dual-mass gyroscope without mechanical linkage to 
oscillate in anti-phase at the same resonant frequency and with the same amplitude. The proposed 
controller also enables the resonant frequencies of the drive and sense axes to be matched precisely by 
continual compensation the frequency deviations, which greatly improves the gyroscope performance. 
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